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ABSTRACT 

We provide an analytical description of the line broadening of H i absorbers in the Lyman-a 
forest resulting from Doppler broadening and Jeans smoothing. We demonstrate that our re¬ 
lation captures the dependence of the line-width on column density for narrow lines in z ~ 3 
mock spectra remarkably well. Broad lines at a given column density arise when the underly¬ 
ing density structure is more complex, and such clustering is not captured by our model. Our 
understanding of the line broadening opens the way to a new method to characterise the ther¬ 
mal state of the intergalactic medium and to determine the sizes of the absorbing structures. 

Key words: cosmology: large scale structure of Universe - quasars: absorption lines - inter¬ 
galactic medium 


1 INTRODUCTION 

The intergalactic medium (IGM) is detected as intervening absorp¬ 
tion in the spectra of background sources (|Gunn&Peterson|1^6^ 
in the form of a ‘forest’ of H i ( |Lynds | ( |1971| l; |Weymann et al. 

I, see e.g. Rauchjj|1998jl^^or a review) and Hen (Jakobsen 
|et al. 1994[ |Davidsen et al. 1996^ absorption lines. This gas 


very highly ionised, x = nni/uH ~ 10~^, by a pervasive back¬ 
ground of ionising radiation produced by galaxies and quasars (e.g. 
|Haardt & Madau|1996f, with hints of an increasingly neutral frac¬ 
tion, X ~ 0.1, above redshift z ~ 7 ( [Mortlock et al. 12011) . At all 
times, it contains the majority of baryons jPukugita et al.|199^ . 

The growth of structure in a cold dark matter (CDM) Universe 
naturally gives rise to a ‘cosmic web’ of voids, sheets and filaments, 
with absorbers arising when the sight-line intersects higher density 
regions, as demonstrated in analytic models JbT (& Davidsen|1997| 
Schaye|[20(^ and simulations ( |Cen et ^|1994[ [Weinberg et al.| 


1996[ Theuns et al.|1998). The temperature of this gas is set by the 


balance between photo-heating and adiabatic expansion and com¬ 
pression, resulting in a well-defined temperature-density relation, 
T = To{p /, where 7 ~ 1 immediately after reionisation, 
and 7— >1-1-1/1.7 asymptotically long after reionisation iHui & 
|Gnedin|1997[|Theuns et al.|1998t . Patchy reionisaton of either H i 
or He II may lead to patchy photo-heating ( [Abel & Haehnelt|1999[ 
|McQuinn et al.|2009| and hence spatial fluctuations in the temper¬ 
ature, but these have not (yet) been detected jTheuns et al.|200^ 
IMcQuinn et al.|20TT) . 

Feedback from star formation is an essential ingredient in 
models of galaxy formation and might impact the IGM ( [Theuns | 
jet al.|20 0T|(. For example, elements synthesised in stars are detected 
in the IGM jCowie et al.|I995t, even at low densities (jSchaye et al.j 
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|2003[ l. If, as is likely, this enrichment is due to galactic winds, then 
the IGM might not be ideal for making cosmological inferences. 
Fortunately, it appears that the effects of such winds are relatively 
small ^Theuns et al. [200^ [McDonald et al.j 2005 |[Viel et al.|2013| (. 
However, current simulations still struggle somewhat to reproduce 
the detected enrichment patterns as a function of density jAguirre| 
[& Schaye[2005[[Wiersma et al.|2010l[2011[ | and hence may under¬ 
estimate the effects of galaxy formation on the IGM. 

Accurate measurements of the evolution of To and 7 are in¬ 
teresting, since they could constrain the epoch of Hi and Heii 
reionisation ( [Theuns et al.[2002) , but also because the formation of 
small galaxies is quenched due to a reduction of atomic cooling H 
|stathiou||1992|l following reionisation, and because photo-heating 
evaporates the gas out of shallow dark matter potentials jOkamot^ 
[et al.|2008) l. This is a crucial ingredient in models of galaxy for¬ 
mation, in particular for dwarf galaxies and satellites ([Benson et al.| 
[2002[ [Sawala et al.|2014k The thermal state of the IGM is also an 
important ‘nuisance’ factor when using the Lyman-a forest to study 
the nature of the dark matter jBoyarskyet^ 2009korto determine 
the amplitude of the UV-background |Boitori^al.[2005[ l. 

Schaye et al. |1999 ' described a method for inferring (To, 7) 


from the line width versus column density (6 — A^hi) scatter plot of 
Lyman-a lines, when the Lyman-a forest is fitted as a sum of Voigt 
profiles for example using VPFITp] This improved and extended the 
work of [Theuns et al.[(r998) who had noted that the small b cut-off 
of the 6-parameter distribution was sensitive to temperature. Schaye[ 
[et ah ([ 2000 ^, and Ricotti et al.[([20Q0|l;[Bryan & Machac"^ ^2000 


and [McDonald et al.[ ( [2001[ l used this to measure the T — p relation 
at redshifts z ~ 2 — 4. 


1 VPFIT is developed by R. Carswell and J. Webb, 
http://www.ast.cam.ac.uk/ffc/vpfit.html 
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Since then other measurements have been performed over a 
larger redshift range, some also based on Voigt profile fitting jRudie| 
|et al.|[20T2), or based on small-scale power spectrum or wavelet 


analysis (Theuns&ZaroubJ2000|JMcDonald^^tay2000yZald^ 


riaga et al.| 2001 [ 

Theuns et al.|2002 Viel et al.|2009 

Lidz et al. 

2010 

IGarzilli et a 

.|20121l, transmission PDF statistics ([ 

lolton et al. 

2008 

IViel et al. 120091 ICalura et al.|2012| 

Garzilli et a 

.| 20121 l and 

more recently on curvature statistics (|Bec 
measurements still have large error bars fc 
even poorer for 7 .|Bolton et al.|(|2008f an 

cer et al .|20 

)r To and cot 
J|Viel et al.| 

"T). Current 
istraints are 
2009J even 


argue that 7 < 1 , which is difficult to understand from a theoretical 
point of view jMcQuinn et al.|2009{|Compostella et al.|2013^ . 

The basis for all these methods is the assumption that line 
widths are relatively directly related to the temperature of the ab¬ 
sorbing gas. However, it is in fact well known that - in the simula¬ 
tions - the absorption lines are invariably wider than just the thermal 
broadening. In addition there is a very large scatter in line widths 
at a given density, in the sense that some lines are much wider than 
the thermal width at a given density. 

In this paper we present a new model for measuring and inter¬ 
preting line widths, which builds on previous work. In our model 
lines are always broader that the thermal width, and we also ex¬ 
plain the origin of the large scatter. We describe how the model can 
be used to infer not just Tq and 7 , but also to constrain the size of 
the absorbers. This paper is organised as follows. In Section 2 we 
begin by contrasting the ‘fluctuating Gunn-Peterson’ description of 
the IGM with one where the absorption is due to a distribution of 
individual absorbers. We do so in order to establish notation and 
explain the motivation for our modelling. We discuss the proper¬ 
ties of these absorbers, and in particular what sets their line widths. 
In Section 3 we introduce the simulations and use them to test the 
basic assumptions of the model. In Section 4 we examine whether 
the method can in principle be applied to observables. We present 
a summary of our conclusions in Section 5. 


2 LYMAN-q scattering in the IGM 

A Lyman-o? photon has a large cross section for exciting a H i atom 
from the ground state to the electronic n = 2 level. When the atom 
falls back to the ground state, the photon is in general not re-emitted 
in the same direction. This scattering of Lyman-a photons out of 
the sight-line to a source is usually called Lyman-a ‘absorption’ 
and is quantified by the optical depth r, such that exp(—r) is the 
resulting transmission - often referred to as the ‘flux’. We begin by 
calculating r for the IGM in two very different approximations. 


2.1 The fluctuating Gunn-Peterson approximation 


The Lyman-a optical depth observed at wavelength A due to a ho¬ 
mogeneous density distribution of neutral hydrogen with proper 
number density nm{z) in a cosmological setting, is given by 


t(A) 



(1 + z)c 
A 



cdz 

A + z)H{z)^ 


( 1 ) 


( |Gunn & Peterson|1965^ , where c is the light speed, a is the Lyman- 
a cross section as function of frequency c/A, Ao ~ 1215.67 A is 
the laboratory wavelength of the H i n = 1 —>■ 2 
H{z) the Hubble parameter at redshift 2 ; (see e.g. 
for an extensive recent review of the relevant physics). Since the 
proper length dr corresponding to a redshift interval dz is dr = 
cdz/ [(1 -f z)H{z)], Eq. ([^ simply states that a column of neutral 


transition, and 
Meiksin|[2009 


hydrogen nmdr contributes a nmdr to the optical depth. Here, Za 
is the redshift of the source. 

Sh'ucture formation results in the growth of the amplitude of 
fluctuations in the density field. This causes the transmission as a 
function of wavelength to be in the form of relatively well defined 
‘absorbers’ l Lynds||1971^, colle ctively referred to as the Lyman- 
a forest ( [Weymann et al.|1981| l. At higher 2 the absorbers blend 
into each other such that regions of near zero transmission are in¬ 
terrupted by small clearings (a Lyman-a jungle), whereas at low 
2 a near unity transmission is occasionally interrupted by a single 
absorber (a Lyman-a savannah). 

When the gas is in photo-ionisation equilibrium, the optical 
depth resulting from a single absorber depends on density and tem¬ 
perature as r oc J.-0.76 y-o.76 (approximately) the 

temperature dependence of the hydrogen recombination rate. Tak¬ 
ing peculiar velocities and thermal broadening into account as well 
results in the ‘fluctuating Gunn-Peterson’ approximation jGnedin| 
|& Hui|19%|[Croft et al.|1998) for the optical depth, 
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where e is the electric charge of an electron, me the electron mass, 
eo is the electric constant, and fiu is the upward oscillator strength 
(usually referred to as the ‘f’ factor). The velocity coordinate v is 
related to the observed frequency through 



(5) 


where tzo = c/Ao is the Lyman a frequency in the lab, 2 is a mean 
redshift of interest. The integral is now over the profile of each indi¬ 
vidual ‘absorber’, the Jacobian \dv'/dx\ takes into account velocity 
gradients, and the sum is over many absorbers. The Gaussian pro¬ 
file expresses (thermal) broadening and neglects the intrinsic line 
profile of the Lyman-a transition. 

The temperature of the IGM is set by the balance between 
photo-ionisation heating of the UV-background, the adiabatic com¬ 
pression, and the adiabatic expansion of the universe. This intro¬ 
duces the temperature-density relation which, for densities around 
the cosmic mean, is approximately a power-law (|Hui & Gnedin| 


T 



To , 


( 6 ) 


where To is the temperature at the mean density (p), and A the 
density contrast. Close to reionisation the gas is nearly isothermal, 
7 ~ 1 , whereas asymptotically long after reionisation, 7 —^ 1 -f 
1/1.7 jHui & Gnedin 1997 Theuns et al. 1998 1 . The 2 < 5 Lyman- 
a forest can be reasonably well described as being a sum of of 
individual absorbers, and we discuss their properties next. 


2.2 Lyman-a absorbers 

The simplest single absorber, a cloud of uniform density, temper¬ 
ature T, and with HI column density Nui, produces a Gaussian 
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absorption line when neglecting the intrinsic Lyman-Q line profile. 


r{v) 


To 
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(7) 


The line centre is at velocity vp — clog(l + z) + Vpec for an 
absorber at redshift z, with peculiar velocity Upec- Equation 0 can 
be obtained from Eq. 0 in the limit that nni \du/dx\~^ has a peak 
at Up with a width in velocity space much less than the thermal 
width b, this is the so called narrow-maximum limit in |Hui et al.| 

iTWf) . 

Several processes contribute to the broadening b of the line. 
The thermal broadening is given by 



with mn the mass of the hydrogen atom and ks Boltzmann’s con¬ 
stant. The second line applies to gas on a temperature-density rela¬ 
tion of the form T — Tp from Eq. 0. 

A more realistic absorber should have at least some density 
structure. Assuming the line centre corresponds to a maximum pp 
in density, a Taylor-expansion in velocity space yields to lowest 
order 


log{p{v)/po) =-{v - vpf/bl. (10) 

Real Lyman-a absorbers may not be strongly peaked, but note that 
what matters is not the density profile of the absorber, but its neutral 
density. Since nni oc for isothermal gas in ionisation equilib¬ 
rium, peaks in neutral gas density are better defined than the max¬ 
ima in density - as we will show using numerical simulations in the 
next section. 

The absorption profile of this more realistic absorber is then 
a convolution of the Gaussians of Eqs. 0 and l |10| l, i.e. another 
Gaussian, with width 

b^=b^T + bl. (11) 

What sets the velocity extent bpl Absorbers at modest density con¬ 
trast are expanding at a rate somewhat smaller than the Hubble 
rate, suggesting that bp ~ XH, with A the physical extent of the 
structure along the line of sight. Based on the papers of |Scha y^ 
( |2001^ and [Gnedin & Hul| ( |1998| l, we expect the physical extent 
of the absorbers to be of the order the Jeans length Aj, defined in 
Eq. ^ below (see also [Miralda-Escude & Rees||1993|( . In the 
model of jSchayej ( [2001^ absorbers are assumed to be locally in 
near-hydrostatic equilibrium, with column density related to den¬ 
sity by 


Nu 


' rim Xj 


( 12 ) 


and hence ‘size’ A,/. Because the dynamical times of low-density 
absorbers are long, the structures are still expanding, and are not 
exactly in hydrostatic equilibrium but are lagging behind. This is 


encapsulated by the ‘filtering scale’ Af = 27r/fcF of Gnedin & Hui 


( [1998^ which depends on the full history of the Jeans length over 
cosmic time and hence depends on the thermal history, not just the 
instantaneous Jeans length. The physical ‘extent’ of the absorber 
in this approximation is then ~ I/Isf, therefore its velocity extent 
is ~ JT {z)/{kF (1 -I- z)) « Aj H{z) /(27r), where Aj is now the 
physical - as opposed to co-moving - extent of the absorber. This 


then motivates us to write 
bp 


r XjH{z) 

JJ —a::— 


(13) 


The factor fj parametrises the time-dependent Jeans smoothing of 
the gas density profiles, and we expect it to be of order unity. 

For gas on a power-law temperature-density, the values of the 
sound speed and Jeans length depend on the density contrast A as 
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(15) 


where p is the mean molecular weight and pt = Pdm + pb is the 
total density, dark matter plus baryons, G is the gravitational con¬ 
stant, Qrn is the matter density. We write the Hubble parameter at 
redshift z as H{z), with Hp = lOOh km s“^ Mpc”^ Rs value to¬ 
day. The numerical values assume p = 0.59, that dark matter and 
gas move together, and for reference we use the jPlanck Collabora-j 
jtion et aT] ( |2013| ) values for the cosmological parameters. Numeri¬ 
cally, the broadening due to Jeans smoothing is thus of order 
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independent of redshift i n the high-redshift appr oximation H{z) oc 
(l + 2)3/3,ThemodelofHui&Rutledgei 1999| uses a very similar 
reasoning to infer line widths; bp is usually referred to as the ‘Hub¬ 
ble broadening’. We demonstrate below that a value of fj ~ 0.88 
describes our numerical simulations well at 2 = 3. 

The total broadening, b, of the absorber is then 


b^ 


bl 

b% 




(17) 


The line-widths of weak absorbers (A<1) are increasingly domi¬ 
nated by Hubble broadening, b^ ~ b^ A~^ oc ^nd hence 

the filtering scales determines their widths. Hubble expansion plays 
a smaller role for stronger lines, with line widths ultimately set by 
thermal broadening alone. Equation jl7| l smoothly interpolates be¬ 
tween filtering and pure thermal broadening. However, we note that 
for typical 2 ~ 3 absorbers with A ~ 3, Hubble broadening still 
plays a significant role, increasing b by almost 70 per cent com¬ 
pared to thermal broadening alone. It might be possible to mea¬ 
sure the extent of absorbers in real space using close quasar pairs 
(see jRorai et al.| ( (2013| l; the impact of Jeans length on quasar pairs 
was also considered in [Peeples et al.|(2010^ ). The extra broaden¬ 
ing discussed here is a longitudinal measure of the Jeans length, as 
opposed to the transverse Jeans length probed with pairs. 

The Gaussian shape of the absorber directly relates column 
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density and width. 




2ro 

(18) 

■j- 'f" _ -7-/2 

t" 

^0 

Am = ^hro, 
o-qc 


(19) 

where r' = dr/dv, t" = cfr/dv^ , and the subscript 0 indicates 
that the quantity is to be evaluated at the (local) maximum of the 


optical depth. We will use these relations below. Equations 
follow directly from Eqs. 0 . 


3 BROADENING OF LYMAN-a ABSORBERS IN 
SIMULATIONS 


3.1 Mock spectra from simulations 


We use the REFERENCE model of the OWLS suite ^Schaye et al.| 
|2010| > to test our model for line broadening. Briefly, the OWLS suite 
is a set of cosmological hydrodynamical simulations performed 
with the gadget-3 code, an improved version of GADGET-2 last 
described by [Spring^ ^2005) , with subgrid models for unresolved 
galaxy formation physics, in particular star formation, and feed¬ 
back by and enrichment from stars, as described in|Schaye & Dalla 


|Vecchia|p008| >, [Dalla Vecchia & Schaye| ( |2008^ and Wiersma et al. 
1 2009^, respectively. Although efficient feedback from star forma¬ 


tion results in outflows from galaxies in this simulation, such winds 
do not have a significant impact on the Lyman-a forest for rea¬ 
sons discussed in |Theuns et al.|j2002| l, see also |Viel et al.| ( [2013| l. 
Most relevant for this paper is that we simulate a A cold dark mat¬ 
ter model using 512^ SPH particles that represent the gas, and an 
equal number of dark matter particles. The computational volume 
is a periodic box of co-moving size 25h~^ Mpc on a side. The 
masses of the gas particles (~ 2 x 10^h~^M q) are sufficiently 
low to obtain numerical convergence for the line widths ( [Theuns | 
jet al.| [T 998) (. Photo-heating of the gas in the presence of an opti¬ 
cally thin UV/X-ray background fromjHaardt & Madau ([2001^ i s 
implemented using interpolation tables from Wiersma et al. 1 2009| l. 

Photo-heating, adiabatic processes, shocks from stmcture for¬ 
mation and feedback from galactic winds and radiative cooling pro¬ 
duce the temperature-density relation shown in Fig. [T^ at redshift 
z = 3. Most of the gas is on a tight temperature-density relation, 
T = To with To « 15000 K and 7 « 1.32 for A<5, and 

this gas is the main reservoir responsible for producing the Lyman- 
a forest. There is gas at significantly higher temperatures at these 
densities, which results from structure formation and feedback. The 
increased importance of radiative cooling decreases the amount of 
hot gas at higher densities, and an equilibrium sequence where ra¬ 
diative cooling and photo-heating balance results in the appearance 
of significant amounts of gas at T ~ 10^ K. This higher density gas 
results in Lyman-limit and damped Lyman-a absorbers (e.g. jAltayj 
|etal.| 2011 ^ . Self-shielding is important for this higher-density gas, 
and the impact of feedback becomes more important ( jAltay et al.j 

[2m3l (. 

Given a snapshot of the simulation at some redshift (we use 
z = Sin this paper), we generate mock Lyman-a spectra using the 
method described injTheuns et al.|( |l998[ ( but using the interpolation 
tables of jWiersma et al. ^2009) to relate total to neutral hydrogen 
density. We generate these spectra at very high resolution using 
pixels of width ~ 2 km s“^, much narrower than any of the ab¬ 
sorption features that appear. Importantly, we neither add noise nor 
instrumental broadening to mimic observed spectra - limiting our 



(a) 



(b) 

Figure 1. Mass-weighted (top panel) and optical depth weighted (bottom 
panel) temperature-density relation in the REFERENCE OWLS run at 2 = 3. 
Values in the bottom panel are evaluated at the locations of maxima in the 
optical depth. Colours are a measure of the density of points in these plots, 
and A = p/(p). In the lower panel the black dots indicate the median tem¬ 
perature in bins of A with eiTor bars computed by bootstrap resampling. 
The black solid line in both panels show the power-law T = ToAT”! 
that best fits the black dots for A < 3. There is a much larger scatter in T 
vs A in the top panel than in the bottom panel, but this shocked gas con¬ 
tributes little to the Lyman-a forest. The lower-envelope of points in the 
T — A relation in the bottom panel mirrors the temperature-density relation 
of the lower-envelope in the top panel, demonstrating that the optical-depth 
weighted temperature at the optical depth maxima of absoiption lines fol¬ 
lows very closely the T — A relation of the photo-heated IGM. 


analysis to comparing our model to idealised observations. Each 
sight line is parallel to a coordinate axis of the simulation volume, 
and we use velocity pixels of size v = H{z)L/[{l + z)N], where 
L = 25h~^ Mpc is the co-moving hox size, and N the number 
of pixels. We calculate t{v), T(v) and p{v), where T and p are 
the optical depth weighted temperature and density, respectively 
(see [Schaye et al.|1999) . Unless explicitly stated, we will always 
use optical depth weighted quantities for variables measured from 
spectra. 

In Fig. we plot values of optical depth weighted tempera¬ 
ture versus density, for pixels that correspond to local maxima in 
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Figure 2. A stretch of a simulated spectrum at redshift 2 : ~ 3. From top to bottom we show the transmission F = (—r), optical depth r, and optical depth 
weighted temperature, T, as function of velocity v. Individual lines are regions of the spectra between two minima in optical depth, they are coloured using 
different colours, with a filled circle at the maximum r for each line. Note that T is not constant over a line. 


the optical depth from these mock spectra. Most points lie close to 
the lower-envelope of the points in the scatter plot, which tracks the 
same T — A relation as the lower envelope in Fig. The scatter 
around this relation is much less in Fig. [Tb because the shocked 
gas that causes the scatter to higher T in Fig. T^contributes little to 
the optical depth. At higher A>10, T decreases with increasing A, 
reflecting the effects of radiative cooling also apparent in Fig.p^ 


3.2 Identifying individual absorbers 

We identify individual absorbers (or Tines’) with regions of the 
spectrum between two minima in the optical depth, dissecting the 
whole spectrum in a set of unique stretches. An example spectrum 
is shown in Fig.[^ with different lines coloured differently, and with 
diamonds at the locations of local maxima in r. Unsurprisingly, the 
temperature varies across a line, and hence assigning a single tem¬ 
perature (or density) to a line is necessarily somewhat ambiguous. 
We could decide to associate to a line the maximum, the (weighted) 
mean or the value corresponding to the maximum of the line for 
any given physical quantity. In the following, where not specified 
otherwise, we will associate to the line the value of the physical 
quantity that corresponds to the optical depth weighted value eval¬ 
uated at the local maximum in optical depth. This method for iden¬ 
tifying lines is not directly applicable to observed spectra in the 
presence of noise. Such spectra should be smoothed to avoid incor¬ 


rectly identifying noise features with small lines, however this may 
lead to missing true weak absorbers. Preliminary analysis shows 
that fitting a Gaussian over a set of contiguous pixels to calculate 
the derivatives that appear in Eq. ID is promising, we will report 
on this in a future paper. 

As demonstrated in Fig. 1, the (real-space) temperature- 
density relation of the IGM in the simulation translates approx¬ 
imately to a power-law (optical depth weighted) temperature- 
density relation for the lines generated in mock spectra. Most lines 
follow r = To oc A^~^ up to A<3, with some scatter to signifi¬ 
cantly higher T, but none scatter significantly below this line. How¬ 
ever for A>3, the median temperature does fall below the power 
law fit: this is because this gas cools radiatively, and this results in a 
decrease in T with higher A for the real-space temperature-density 
relation as well. Note further that the T — A relation is only ap¬ 
proximately a power-law, with evidence for some curvature even 
over a relatively narrow A range. Finally, note that in the absence 
of noise, we see that the majority of lines have 0.1<A<10. 

The optical depth profile of several of the lines in Fig. 0 ap¬ 
pear reasonably Gaussian in shape near their maxima, and we will 
investigate to what extent their widths follow Eq. However it 
is also clear that some lines blend with nearby lines, and we will 
show below that this ‘clustering’ impacts how well we can infer 
the underlying T — A relation from spectra - even in the absence 
of noise. 
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Figure 3. Velocity extent of absorbers in real space versus density 
contrast A for lines identified in mock spectra without peculiar velocities. 
Colours indicate density of points in the plot. The 5th, 10th, 90th and 95th 
percentiles of in bins of A are shown as symbols with bootstrap er¬ 
rors. Line widths decrease with increasing density, opposite to the case of 
thermal broadening. This behaviour is captured by the solid black line, the 
broadening bp from Eq. dn for/j = 0.88, which fits the lower percentiles 
of the b/^ — A relation well in shape and amplitude for A^3. 



Figure 4. Total broadening b^ = —2 tIt" of absorbers at line centres as a 
function of density contrast A for lines identified in mock spectra. Colours 
indicate density of points in the plot. Black dots with bootstrap error bars 
denote the location of the 10th lower percentile of b in bins of A. The dotted 
red line is the thermal width from Eq. I^for the temperature-density relation 
from Fig.^j lines are invariably wider than bj^ except at A^IO. The solid 
black line is the model of Eq. {3 with fj = 0.88, for the same T — A 
power-law relation. It describes the black dots very well up to A ~ 3, 
above which the simulation’s T — A relation is not well fit by a power- 
law. The dashed magenta line is again the model of Eq. (3 but now using 
the measured T — A relation rather than a power-law fit: it captures the 
down-turn in b for A>3 due to the onset of radiative cooling. 


We suggested that line broadening has a contribution bp due 
to the spatial extent of the absorber as described by Eq. ( |13^ . To 
test this directly, we should plot bp versus A. However, maxima 
in T do not correspond directly to maxima in A because of pe¬ 
culiar velocities. Therefore, to unambiguously identify a line with 
a density structure, we generate mock spectra setting peculiar ve- 


Figure 5. Histogram of line broadening = —2t/t" for simulations with 
and without peculiar velocities, shown as solid black and dashed red lines, 
respectively. The 10th percentiles of narrowest lines are indicated by verti¬ 
cal lines. The impact of peculiar velocities on the distribution of b values is 
small. 


locities to zero. We later show that this has very little effect on 
the line broadening (see also Fig. |^. For each line identified in 
these mock spectra, we can now unambiguously identify the cor¬ 
responding density stretch that gives rise to that line. However, the 
shape of the density structure of lines is often complex and it is 
not immediately apparent what extent to associate with a given 
line. We decided to use the following prescription. For each line 
we determine its start and end velocity, vi and V 2 , the location 
(Vmax) and height (Amax) of the maximum density, and the integral 
I = ^{v) dv. We now determine the value of 6 a of the Gaus¬ 


sian profile, G(A) = Amax exp{—{v — Vmax)^/6A). for which 
G{v) dv = I. In other words, we define the extent of the ab¬ 
sorber, 6 a, as the width of the Gaussian that has the same integral 
and maximum value as the line itself. If the line were a Gaussian, 
then 6a would simply be its width. Further details on the fitting 
procedure can be found in Appendix [B] We expect that 6a ~ bp 
from Eq. 


In Fig 


Gil- 

PigH' 


Jwe plot 6^ versus A for all lines identified in mock 
spectra generated while ignoring peculiar velocities. As before, 
colour encodes the number density of lines in this plot, and 5th, 
10th, 50th and 95th percentiles are plotted as symbols with boot¬ 
strap error bars. Values of 6 a increase with decreasing A - the op¬ 
posite of what is expected from thermal broadening. The solid line 
indicates 6^ from Eq. 1 13 i with fj — 0.88 applied to gas with the 
T — A relation T = 15000 K which, as we showed earlier, 

fits the actual temperature-density relation of the simulation well 
for 0.1<A<3. Fig.|^shows that both the 5th and 10th percentiles 
follow the solid line well over the same A range. This implies that 
bp, with fj = 0.88, is indeed a good estimate of the width of the 
absorbers in real space. 

Given that our estimate of the broadening of the density pro¬ 
file works well, we plot the total line broadening (including peculiar 
velocities) versus density contrast, and compare it to the model of 
Eq. (Fig.g. As before, the colour represents the density of 
points in the plot. We compute 6 from the curvature of the optical 
depth profile t{v) using Eq. 1 18 ; fitting a Gaussian to the pixels 
close to the maximum yields very similar values of 6. Values of 
6 — A for individual lines show a larger scatter in 6 at given A, but 
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Figure 6. Optical depth r {top panel) and corresponding density contrast A {bottom panel) for mock spectra for the simulation with the original temperature- 
density relation {solid red line) and for the case where T = 100 K ever ywhere (dotted black line). Clustering of density peaks in real space makes some 
of the absorption lines much wider than the estimate of b from Eq. liz}’ for example the wider lines around v ~ 800 km s ^ and v ~ 1000 km s 
Thermal broadening and Hubble broadening together smooth the individual peaks seen in the underlying density field into one wide line, whose width cannot 
be described accurately by eq. [T7] 


there is a well-defined lower envelope below which few lines ap¬ 
pear. This lower envelope is traced by the 10th percentile of b values 
in narrow bins of A, plotted using black symbols with bootstrap 
errors. The red dashed line indicates the thermal broadening, br, 
from Eq. using the power-law fit to the temperature-density re¬ 
lation measured in the simulation. For A<10, the measured broad¬ 
ening along the cut-off is significantly wider than br, with the dis¬ 
crepancy between the two increasing with decreasing A. The solid 


black line shows -|- 6p as a function of A from Eq. 117 again 
using the best-fitting power-law temperature-density relation, and 
adopting the value fj — 0.88 found above. It captures accurately 
the dependence of the broadening on A for the lower envelope of 
the absorption lines traced by the 10th percentile for A<3, in par¬ 
ticular reproducing the measured upturn in b^ for A<1. Because 
the solid black line assumes that the relation between temperature 
and density is a power law, it does not describe the simulated b — A 
relation well above A ~ 3. The dashed black line uses the me¬ 
dian value of the temperature at a given density for absorption lines 
when calculating the Jeans smoothing term bp in Eq.^l3[l, rather 
than a power-law. It does much better in capturing the downturn in 
6 —A at A>10. We conclude from this plot that the model in which 
the line broadening is a combination of thermal and Jeans smooth¬ 
ing, Eq. describes the lower-envelope in the b — A plane well. 
Note that deviations of the true temperature-density relation from 


a power-law, and Jeans smoothing, cause the T — A relation to 
be non-monotonic, with a minimum value of b for A ~ 1 in this 
particular simulation at this particular redshift. 

Peculiar velocities are not the cause of the appearance of lines 
that are much wider than the broadening b computed from Eq. jl7| l. 
We compare the distribution of b values for lines identified in mock 
spectra with and without peculiar velocities in Fig.[^ The peculiar 
velocities do not have a large effect on the broadening, if anything 
they make lines slightly narrower, see also |Theuns et al.| ( |2000| l. 

Our analytic expression for broadening, Eq. GZ). provides 
only a lower limit on the width of a collection of lines in the Ly¬ 
man a forest. The additional broadening can be partially attributed 
to the underlying density structure in the absorber jHui & Rutledge| 
|1999| l, for example in terms of the angle under which the sightline 
intersects the filament that corresponds to the absorber. However, 
another important contribution that is not easily captured by linear 
theory is the clustering of density peaks. We illustrate this in Fig.|^ 
comparing a mock spectrum with the original T — A relation with 
the same spectrum now computed assuming T = 100 K every¬ 
where. We see that some individual absorption lines identified in 
the solid red spectrum have several underlying density maxima. If 
the temperature of the gas is low, as in the case of the black dotted 
spectrum, then these individual density maxima give rise to their 
own absorption lines. But if the gas is sufficiently hot, then thermal 
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Figure 7. The neutral hydrogen column density of absorption lines versus 
density contrast for the REFERENCE model at 2 = 3; colours are a measure 
of the density of points in the plot. Black dots denote the 10th percentiles of 
column density in logarithmic bins of A. The solid black line corresponds to 
the model of Eq. {ig. using the value of T appropriate for the REFERENCE 
model with values of Tq, 7 and fj taken from Fig.[^ and for /jv = 0.3. 
This value of is the result of best-fitting Eq. dD to the 10 th percentiles, 
for A 0.3. This provides an excellent fit to the results from the simulation 
over 4 orders of magnitude in column density. 


and Hubble broadening can merge multiple individual components 
into one much wider line. The width of that absorber is not de¬ 
scribed well by our model. We conclude that it is the clustering of 
density maxima that gives rise to the lines that are much wider than 
expression b from Eq. (TtJ. Another way to demonstrate this is to 
count the number of individual density peaks for each absorber: we 
find that absorbers with widths close to that from Eq. jl7| l typically 
have only one density maximum, whereas those that are signifi¬ 
cantly wider correspond to several underlying peaks. 


4 INFERRING THE TEMPERATURE-DENSITY 
RELATION FROM OBSERVABLES 


Our description of line broadening so far used b and A to relate 
temperature to density, yet density contrast A is not directly mea¬ 
surable. Here we investigate whether we can use the column density 
Nui of the line instead. Following the reasoning of |Schay^ ( |2001| l, 
we make the Ansatz that the column density, density, and extent of 
the line are related by 

Am ~ nni/jAj . (20) 


For highly-ionised gas in photo-ionisation equilibrium 
n-HI 

- = aniinei 

riH 


X 10 “® (. 


Vl.2- 10 

-0.76 / 


riH 


' cm~^ 
\ -1 


r) 


where 

j, \ -0.76 

Wk) 

ao = 4 X 10”^^ cm^s”^ , 


anil 


ao 


( 21 ) 


( 22 ) 


is the value of the case-A recombination coefficient as a function of 
temperature, rie is the electron number density, nn is the total hy- 


A 



10 '^ 10 ^^ 10 ^" 10 '= 10 '® 10 " 

AHi[cm^ ] 



10th percentiles in space 



Figure 8. The line-width column density — A^hi relation in the REF¬ 
ERENCE OWLS model at redshift 2 = 3, with broadening determined from 
Eq.dD and column density from Eq. G3- The top axis gives the value of 
the density contrast that corresponds to Nm from Eq. (D- Colours denote 
the density of points in this plot, black dots denote the 10 th percentile of 
in bins of column density with bootstrap error bar s. Th e dotted red line 
is the thermal broadening b^ from Eq. using Eq. l |23| to relate column 
density to density contrast. The solid bla^ line is our model for broadening 
of Eq. j25t using the values To and 7 from Fig[T^ fj = 0.88 from Fig.||] 
and /tv = 0.3 from Fig.[7| The black line connecting magenta dashes is 
Eq. j25t using the T-A relation measured from the simulation, rather than 
a power-law fit. The analytic formula describes the minimum line width at a 
given column density accurately, and works even better if the actual T — A 
relation is used rather than a power-law fit. 


drogen number density, and F the photo-ionisation rate. Combining 
these with the expression of the Jeans length from Eq. yields 


Nm = 


3^ 


GTtt 
ks lO'K 


1/2 


Gj ao 


Hi 


fN fj f T 


0.3 0.88 (,104k ) 


T 

ICHK 
-0.26 


r (GmH)2 

-0.26 


A3/2 


(23) 


where numerically 
No = 1.4 X 10 '® cm 

Hi, " " 


0.04825 


H, 


10-12 s-i 
- 1/2 


1 + z 


9/2 


0.307 


h 


0.6777 


(24) 


Here, fjv is the proportionality factor implicit in Eq. l |20| l which we 
show below is ~ 0.3, and the numerical value for No is for very 
highly ionised (fi = 0.59) primordial gas with the [Planck Collabo^ ] 
[ration et ar] ( |2013) cosmological parameters, and Y is the primor¬ 
dial helium abundance. Eor a power-law temperature-density rela¬ 
tion, T — To A'*'-!, the hydrogen column density depends only 
weakly on temperature, Nm oc and it scales with density 

contrast as /\it 6 -o .267 compare the results from Eq. (23 i to 
our simulation in Fig. [7] where we find the best fit of the model to 
the 10th percentile distribution of the Nm versus A, the justifica¬ 
tion of this choice is in Appendix]^ The model fits very accurately 
the results from the simulation over nearly four order of magni¬ 
tude in column density, as also shown in |Rahmati et al.[ ( [2013^ and 
[Tepper-Garcfa et ar]j2012| . 
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Finally, we combine Eq. l |23| l with Eq. l |17| > to obtain a relation 
between the two observables quantities, namely the column density 
A^hi and the lower-envelope of line widths at a given column den¬ 
sity, b, in terms of underlying temperature-density relation of the 
gas, as 


bx 

b^ 

N 


'^ksTp ^{ 7 _ 1 )/{ 1 , 76 - 0 . 267 ) 

m 


1 + 0.75 


f Ar-l/(l'76-0.267) 

\0.88j 


NuA f To \ 
No ) 




(25) 


This relation is compared to the b — Nm distribution of the simula¬ 
tion in Fig.[^ Our model for line broadening does well in predicting 
the minimum line width as a function of column density, in partic¬ 
ular predicting correctly that line width is not a monotonic function 
of column density. At low A, line widths increase as Jeans smooth¬ 
ing becomes more important, whereas at high A, line widths de¬ 
crease due to radiative cooling. The underlying T — A relation of 
the simulation is only approximately a power law: using the rela¬ 
tion from the simulation directly in Eq. l |23^ improves the agree¬ 
ment between model and simulation further, in particular it might 
be possible to measure the temperature of denser gas, A ~ 10. We 
illustrate the sensitivity of the b — Nui relation to small changes in 
To and 7 , as well as the fitting parameters fj and /jv, in Appendix 
C. 

Our model overestimates line broadening of very weak lines, 
cm~^, yet it correctly predicts line widths at low den¬ 
sity contrast (see Fig. |^. In Appendix A we show that this is be¬ 
cause line width biases how well a line of given column can be 
associated with a given density contrast. In practise these lines are 
too weak to be detected in observed spectra. 


5 CONCLUSIONS 

We have presented an analytic description of the broadening of 
lines in the Lyman-a forest, in terms of the temperature, Jeans 
smoothing and density of the absorbers. Identifying individual ab¬ 
sorbers in a spectrum as spectral stretches between two minima in 
the optical depth, we calculated the line width, 6 , from the sec¬ 
ond derivative of the optical depth at the local maximum in r, 
= - 2 r/r" (Eq|l 8 |. We argued that the thermal width of a 
line, br, combined with its width due to Jeans smoothing, bp, sets a 
lower-limit on the total line width, b^>b^ + b^. We derived an ex¬ 
pression for bp that depends on the temperature T and the density 
contrast A of the absorber, which smoothly interpolates between 
the filtering length interpretation of line-widths from |FIui & Gnedin| 
(|1^97b at low A, and the near hydrostatic equilibrium model of 
|Schaye| ( [2001^ at higher A. This model gives a good description of 
the lower values of b at given A (Fig.|^ . We also discussed the ori¬ 
gin of the lines that are significantly wider than the lower envelope, 
arguing that they result from the clustering of several individual 
density maxima that cannot be distinguished in the spectrum due 
to thermal broadening and Jeans smoothing. Such clustering is not 
captured by our quasi-linear model. 

Based on the model of |Schaye] ( |2001^ , we write an expres¬ 
sion relating column density to density contrast. Combining this 
with our model for line broadening yields an analytic expression 
for the lower envelope of 6 as a function of column density, Nm 


(Eq.|25|>. This model describes the lower envelope much better than 
one based on thermal broadening alone (Fig.[^ . 

In our interpretation of b, line widths of weaker lines are dom¬ 
inated by the filtering length, which in principle depends on the 
thermal history of the IGM and not just on its instantaneous state. If 
this is true, then not taking this into account may affect the interpre¬ 
tation of the b—Nm cut-off in terms of the underlying temperature- 
density relation of the gas. 

Our analytic description of the line broadening opens the way 
to a new method for determining the IGM thermal state, and in 
particular for measuring the ‘longitudinal’ size of the absorbers (the 
‘filtering length’), that is complementary to what can be measured 
from close QSO pairs. This is particularly interesting, because for 
the first time we have a method to estimate the sizes of the clouds 
from a single line of sight. 
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APPENDIX A: BIASES IN THE Ahi-A RELATION 

In Fig. |A1| we plot A-Ahi for lines, with Ahi computed from 
Eq. l|19|), and individual lines coloured with the value of the broad¬ 
ening b, calculated from Eq. l |18| >. Line width biases the value of 
column density, in the sense that the analytic model of Eq. 123 1 
(black line in the figure) yields lower values of Ahi than Eq. 119 1 
for broad lines. The effect of this bias on the Nm-b relation is il¬ 
lustrated in Fig. |A2| Our model for line broadening (black line) 
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Figure Al. Density contrast - column density A-A^hi relation for lines 
from the REFERENCE model, coloured according to their line width, 6; black 
dots denote the 10th percentile lowest value of column-density at given den¬ 
sity contrast. Line width b and column density A^hi calculated from 
Eqs. jl8|19) , respectively. The solid black line is the analytic relation from 
Eq. j23t - At given A, broader lines have higher column density, illustrating 
that line width biases the column density associated with a given density 
contrast. 


A 





Figure A2. Column-density line width N^i-b relation for the REFERENCE 
model, colours are a measure of the density of points in the plot; the black 
line is the model illustrated in Fig.[^ the doited red line is the thermal width. 
Black dots denote the 10th percentile lowest values of b for bins in column 
density, magenta triangles denote the 10th percentile lowest values of b for 
bins in density contrast A, using Eq. to convert A to Nm. The bias 
that line width introduces in this relation causes the black dots and red dots 
to diverge at low A^hi ■ 


works better when lines are characterised by their density contrast 
(red dots) than by column density (black dots), especially for weak 
lines. This is because relation Eq. 1 23 i is fitted to narrower lines. 

Here we also show that Eq. ' 19' is equivalent to computing 
Vhi with the integral of the optical depth, as it is often done in the 
literature. In Fig.|A3|we show A-A^hi, with Nm computed from 


1 

Nm = - / rdv, (Al) 

^oc 


where ui and V 2 are the extremes of the line. If we compare Fig.|A3| 


Figure A3. Same as Fig. 1^ but with Nm computed from Eq. jAl) . For 
the low-density lines the neutral hydrogen column density deviates from 
Eq. due to the blending of the faint lines. 



A 


Figure A4. Same as Fig. |A1| and |A3| but with Afjji computed from 
Eq. jA2) . The distribution A-Ahi looks similar to the one shown in 

Fig.[^ 


with Eig. |A1[ we can see that Eq. | |A1^ underestimates the neutral 
hydrogen column density for the low-density lines. 

When lines are blended, the column density computed from 
Eq. ^ may become smaller and smaller as a shorter velocity 
interval is associated with the line (i.e. when vi ~ V 2 ). We can 
estimate the contribution of these wings, and compute a corrected 
column-density 

1 /■°° 

Nm = - / rdv, (A2) 

o-QC J_^ 

by extrapolating a Gaussian fit to the line. Such an extrapolated 
density is closer to what is done in Voigt-profile fitting. Using 
Eq. l |A2| > results in a 6 — Ahi distribution which is very similar 
from that obtained in Fig.|Al| as shown in Fig.|A4| 
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Figure Bl. Density profile of absorbers in the absence of peculiar velocities {black line) for a mock spectrum calculated from the REFERENCE model; dashed 
blue line is the scaled optical depth for comparison. Red lines are Gaussian profiles from Eq. where the width is such that the integral of the density 
profile of the Gaussian is the same as the one of the multi-peaked absorber. The widths of the red Gaussians are often a good approximation to the extent of an 
absorber. 


APPENDIX B: DETERMINATION OE fj 

In this Appendix we illustrate that our estimate for the extent of ab¬ 
sorbers in real space indeed is a good approximation for their actual 
size. We begin by neglecting peculiar velocities (i.e. generate mock 
spectra after zeroing all peculiar velocities), so that we can unam¬ 
biguously identify the real-space density structure that gives rise to 
a given line; we demonstrated in the main text that peculiar veloc¬ 
ities have little impact on line-widths. Often individual absorption 
lines correspond to more than one density peak. Such clustering of 
peaks sometimes gives rise to very wide lines, as we argued in the 
text. 

In Fig. BTIwe show how we have computed the broadening 
6 a in Fig.lS from the density profile of the multi-peaked absorber. 
The red lines are Gaussian profiles as in Eq. with Vo and po 
respectively the velocity and the density of the local maximum in 
optical depth. The width 6 a is such that the integral of the Gaussian 
density profile is the same as the integral of the density profile of the 
absorbers. Even though there is often no unique way of associating 
a ‘size’ with such complex absorbers, we believe 6a is often a good 
approximation of the extent over which the density contrast of an 
absorber is significant. 


APPENDIX C: DEGENERACY OE PARAMETERS IN THE 
LINE BROADENING RELATION 

For completeness we illustrate the sensitivity of Eq. GZl to changes 
in the T — A relation (varying Tq and 7 ; Fig[CT]l, and varying the 
fitting parameters fj and /jv (Fig.|C2^ 


A 



10 ^^ 10 ^^ 10 ^" 10 ^^ 10 ^® 10 ^^ 

AHi[cm^ ] 



^ 10th percentiles in space 



Figure Cl. Same as figure]^ The dashed (dashed-dot) lines are the model 
from Eq. |25| with To ( 7 ) changed by ±20%. 


A 



10 ^^ 10 ^^ 10 ^'* 10 ^® 10 ^® 10 ^^ 

AHi[cm^ ] 



^ 10th percentiles in space 



Figure C2. Same as figure]^ The dashed (dashed-dot) lines ai‘e the model 
from Eq. |25| with fj {Jn) changed by ±10%. 

































































